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ABSTRACT 

This  research  analyzes  the  complex  permittivity  and  permeability  of  thin  samples 
of  materials  with  large  air  gaps  between  the  material  sample  and  outer  conductor  of  a 
coaxial  waveguide.  Assuming  axially  symmetric  air  gaps,  higher  order  fields  are  excited 
in  the  sample  region  which  must  be  accounted  for  due  to  the  power  coupling  into  these 
higher  order  fields.  This  work  builds  upon  previous  research  which  analyzed  the  effects 
of  smaller  air  gaps  on  the  extracted  material  parameters  of  both  dielectric  and  magnetic 
material  samples. 

This  research  utilizes  the  modal  method  developed  in  the  previous  research  to 
calculate  theoretical  scattering  parameters  of  the  sample  region  as  a  function  of 
permittivity,  permeability,  and  frequency.  Complex  root  search  algorithms  (either  a 
Newton  or  Eevenberg-Marquardt)  iterate  the  values  of  permittivity  and  permeability  until 
the  difference  between  the  theoretical  and  measured  scattering  parameters  is  minimized. 
The  final  theoretical  permittivity  and  permeability  values  are  determined  the  be  the  actual 
properties  of  the  material  sample. 

Accurate  results  were  obtained  for  a  dielectric  sample  which  filled  only  22%  of 
the  radial  distance  between  inner  and  outer  conductors  when  considering  20  modes.  It  is 
shown  that  the  solution  appears  to  converge  after  consideration  of  only  10  modes. 

When  utilized  with  20  mil  samples  of  the  magnetic  material  WaveX  filling  only 
12.5%  of  the  distance  between  inner  and  outer  conductors,  the  modal  method  returned 
very  poor  results  for  permittivity  and  relatively  accurate  permeability  results.  The 
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WaveX  samples  were  cut  from  flat  sheets  of  material  and  “wrapped”  around  the  inner 
conductor,  likely  resulting  in  errors  which  are  not  accounted  for  in  the  modal  method 
development.  Additionally,  the  material  exhibited  strong  anisotropic  tendencies  which 
also  were  not  accounted  for  in  the  modal  method  development.  It  appears  that  the  higher 
order  modes  evanesce  much  more  rapidly  than  in  the  dielectric  material  case,  resulting  in 
less  impact  to  the  system. 

Two  sources  of  error  likely  significantly  degraded  the  effectiveness  of  the  modal 
method  for  the  thin  WaveX  samples.  The  first  source  is  measurement  error  which  is 
attributed  to  the  discontinuities  created  when  wrapping  the  material  around  the  conductor. 
This  likely  invalidated  the  assumption  of  axial  symmetry  which  would  lead  to  a  more 
complex  field  development.  The  second  error  source  is  the  assumption  of  a  simple  media 
(i.e.  linear,  isotropic,  and  homogeneous).  It  was  discovered  by  the  end  of  this  research 
that  the  WaveX  is  in  fact  highly  anisotropic.  This  fact  would  lead  to  an  entirely  different 
field  development  which  would  greatly  impact  the  results  of  the  modal  method. 
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LOW  FREQUENCY  MATERIAL  CHARACTERIZATION  OF  THIN 


SUBSTRATES  IN  A  COAXIAL  TRANSMISSION  LINE 

/.  INTRODUCTION 

All  materials  contain  certain  intrinsic  properties  which  determine  their  response  to 
electromagnetic  fields  in  the  environment.  In  today’s  high  technology  environment, 
electromagnetic  fields  are  ever  present  from  applications  ranging  from  personal 
communications  to  defense  applications  such  as  threat  radars.  This  proliferation  of 
electromagnetic  energy  drives  the  need  to  continue  developing  new  materials  with 
specific  responses  to  these  electromagnetic  fields.  The  electromagnetic  characteristics  of 
these  materials,  primarily  permittivity  and  permeability,  are  the  core  to  understanding  the 
material’s  response  to  an  induced  electric  and  magnetic  field  respectively. 

A  classic  example  of  an  application  in  which  this  knowledge  is  critical  is  the  field 
of  electromagnetic  interference  and  compatibility  (EMI/EMC).  Most  modem  electronic 
devices  work  in  close  proximity  to  each  other  and  therefore  can  be  susceptible  to 
electromagnetic  interference  from  each  other.  To  overcome  this  interference,  certain 
materials  may  be  incorporated  into  the  device  to  act  as  a  shield  to  prevent  unwanted 
external  energy  from  penetrating  the  device.  EMI/EMC  can  prove  to  be  a  very  large 
threat  to  complex  systems  such  as  aircraft  due  to  their  operating  environment  and 
onboard  high  power  communications  and  navigation  equipment.  Today’s  military 
aircraft  employ  a  variety  of  materials  in  their  design  which  serve  as  shields  to  protect 
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vital  systems  such  as  avionics  and  munitions  from  this  external  electromagnetic 
interference.  These  materials  are  designed  to  absorb  some  of  this  energy,  preventing  it 
from  reaching  the  critical  systems,  and  lower  the  amount  that  might  be  scattered  back 
into  the  environment. 

These  materials  are  becoming  increasingly  thin  and  flexible  thanks  to  furthered 
understanding  in  materials  science  and  advances  in  manufacturing  technologies.  One 
problem  created  with  these  more  complex  materials  is  the  performance  of  traditional 
methods  used  to  characterize  the  permittivity  and  permeability  of  the  material.  This 
drives  the  need  to  develop  novel  techniques  that  quickly  and  accurately  determine  these 
properties  in  a  non-destructive  manner. 

1.1  Problem  Statement 

The  sponsor  of  this  research,  the  Air  Force  Research  Laboratory,  Sensors 
Directorate  (AFRL/RY),  expressed  interest  in  the  electromagnetic  characterization  of  thin 
materials  in  low  frequency  environments.  Many  techniques  such  as  transmission  lines, 
cavities,  free-space  systems  and  striplines  currently  exist  to  measure  a  material’s 
electromagnetic  properties  at  varying  frequency  ranges.  Low  frequency  characterization 
carries  inherent  challenges,  most  notably  the  sample  size  required  for  accurate 
measurements  at  these  frequencies.  Many  of  these  techniques  previously  mentioned 
require  impractically  large  sample  sizes;  therefore,  many  traditional  methods  are 
eliminated  from  consideration. 

One  technique  which  shows  promise  and  is  the  focus  of  investigation  in  this 
research  is  the  TEM  coaxial  transmission  line/waveguide.  The  design  of  the  TEM  line 
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allows  for  low  frequency  interrogation  using  relatively  small  sample  sizes.  Additionally, 
the  coaxial  line  has  a  normal  (i.e.  radial)  or  vertically  polarized  electric  field  to  samples 
which  are  wrapped  around  the  inner  or  outer  conductors.  This  allows  for  a  strong 
traveling-wave  field  interrogation  of  a  thin  material  wrapped  around  either  conductor. 
Additionally,  the  coaxial  line  provides  an  axially  symmetric  environment  which  reduces 
the  complexity  of  the  fields  present  in  the  guide. 

Most  waveguide  measurements  use  the  Transmission/Reflection  (TR)  method  to 
extract  the  material  parameters  from  the  sample.  This  method  utilizes  a  Network 
Analyzer  (NWA)  to  measure  the  scattering  parameters  of  the  material  contained  in  the 
waveguide.  An  algorithm  such  as  Nicolson-Ross-Weir  (NRW)  is  widely  used  by  the  TR 
method  as  a  closed  form  method  to  extract  the  material  properties  from  the  measured 
scattering  parameters  [15-16].  The  fundamental  assumption  of  the  NRW  algorithm  is 
that  the  material  sample  fully  fills  the  cross  section  of  the  guide  which  restricts  the  fields 
in  the  system  to  a  single  dominant  mode. 

Typically,  samples  are  machined  which  radially  fill  the  line  and  have  a  relatively 
small  thickness  or  length.  At  low  frequencies,  the  wavelengths  of  electromagnetic  waves 
can  become  very  large.  For  example,  the  wavelength  of  an  electromagnetic  wave  at  100 
MHz  is  approximately  3  meters.  Typically,  most  samples  are  machined  to  be  no  more 
than  a  few  inches  thick;  therefore,  they  appear  very  “electrically  thin”  to  low  frequency 
incident  waves.  One  way  to  practically  lengthen  the  sample  to  accommodate  for  this 
electrical  thinness  is  to  manufacture  thicker  but  less  radially  wide  samples.  This  longer 
sample  allows  for  more  of  the  low  frequency  wave  to  impact  the  sample  theoretically 
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improving  the  accuracy  returned  by  the  scattering  parameters.  One  major  downside  of 
this  approach  is  the  introduction  of  radial  air  gaps  between  the  inner  or  outer  conductors. 

Air  gaps,  such  as  those  in  the  setup  shown  in  Figure  1.1,  cause  the  excitation  of 
higher  order  modes  which  invalidates  the  single  dominant  mode  assumption  made  by  the 
NRW  algorithm.  Energy  from  the  incident  mode  is  coupled  into  these  higher  order 
modes  at  the  material  interface,  most  of  which  are  highly  evanescent.  These  higher  order 
modes  are  scattered  back  and  transmitted  through  the  material  to  the  opposite  side.  The 
detection  ports  of  the  NWA  measure  only  the  dominant  mode;  therefore,  the  energy 
coupled  into  the  higher  order  modes  is  not  accounted  for  in  the  reported  scattering 
parameters.  In  this  case,  a  closed  form  solution  such  as  the  NRW  algorithm  may  result  in 
dramatic  errors  in  the  extracted  properties  due  to  the  unaccounted  energy  resident  in  the 
higher  order  modes. 

This  primary  approach  considered  by  this  thesis  is  to  perform  a  modal  analysis  of 
the  system  to  produce  theoretical  scattering  parameters  which  account  for  the  energy 
coupled  into  the  higher  order  modes.  These  theoretical  values  are  used  along  with  the 
measured  values  to  determine  the  actual  permittivity  and  permeability  values  of  the 
material.  Other  approximate  methods  such  as  a  perturbation  analysis  or  closed  form 
correction  factors  have  shown  promise  in  previous  air  gap  research  and  may  also  be 
worthy  avenues  of  pursuit.  The  details  of  the  modal  method  and  other  possible  methods 
are  discussed  in  detail  later  in  Chapter  3. 
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1.2  Limitations 


The  TR  method  is  very  useful  in  determining  the  complex  permittivity  and 
permeability  of  a  material.  For  each  unknown,  one  independent  measurement  is  required; 
therefore,  to  determine  both  permittivity  and  permeability,  two  independent 
measurements  are  required.  Typically,  the  two  independent  measurements  used  are  the 
Reflection  (5'jj)  and  Transmission  (5'2i)  coefficients.  In  the  system  shown  in  Figure  1.1, 
the  thin  material  fills  only  a  very  small  portion  of  the  cross  sectional  area  of  the  coaxial 
line.  In  this  case,  a  strong  reflection  response  from  the  material  is  not  expected  due  to  the 
small  cross  section  of  material  to  perturb  the  field.  In  the  event  that  the  reflection 
response  is  inadequate  to  accurately  perform  the  modal  analysis  technique,  an  additional 
independent  S2]  measurement  is  needed  to  provide  the  necessary  equations  to  determine 
both  parameters  of  interest.  An  independent  5*2 1  measurement  may  be  obtained  by 
measuring  samples  of  varying  lengths  or  radial  thickness. 

An  additional  limitation  of  this  research  is  the  availability  of  accurate,  well 
quantified,  permittivity  values  for  the  WaveX  material  used  in  the  measurements.  Some 
WaveX  measurements  have  been  conducted  by  the  Georgia  Tech  Research  Institute 
(GTRI)  and  provide  permittivity  and  permeability  data  for  varying  frequency  ranges. 
Neither  set  of  measurements  are  statistically  well  quantified,  as  the  reported  properties 
for  different  samples  may  vary  extensively.  It  is  anticipated  that  the  results  from  the 
modal  method  utilized  in  this  research  will  coincide  with  the  measured  values  from  GTRI 
shown  in  Appendix  D  within  a  reasonable  margin  of  difference. 
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1.3  Scope 

Other  methods,  sueh  as  mierostrip  probes  eould  be  used  to  approaeh  the  problem 
of  low  frequeney  measurements  of  thin  materials.  Previous  researeh  developed  the  Full 
Wave  Modal  Method  for  a  eoaxial  transmission  line  with  inner  or  outer  air  gaps  [8].  This 
thesis  builds  upon  this  method  and  its  effeetiveness  in  extraeting  the  eomplex  permittivity 
and/or  permeability  of  very  thin  dieleetric  and  magnetie  materials  from  the  measured 
scattering  parameters.  The  method  is  verified  first  using  samples  of  a  low  loss  dielectric 
material  and  then  expanded  to  a  lossy  magnetic  material.  The  magnetic  material  (WaveX 
WX-A020)  is  20  mils  thick  and  is  specifically  designed  for  EMI  absorption  at  low 
frequencies.  Considering  the  possible  limitation  of  the  S'jj  measurement  described  in 
Section  1 .2,  multiple  length  samples  were  used  to  obtain  independent  measurements. 

The  lower  bound  of  the  frequency  range  is  set  to  50  MHz  which  is  slightly  above  the  10 
MHz  lower  bound  of  the  Network  Analyzer.  In  theory,  the  modal  method  should  prove 
valid  for  all  frequency  ranges  at  which  measurements  can  be  conducted.  Other  methods 
such  as  striplines,  perturbation  techniques,  and  air  gap  corrections  may  be  investigated 
for  purposes  of  comparison  to  the  coaxial  line  modal  method. 
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Figure  1.1:  Low  Frequency  Coaxial  Transmission  Line  Concept 
1.4  Thesis  Organization 

Chapter  2  contains  an  overview  of  electromagnetic  fundamentals  starting  with 
Maxwell’s  equations  and  moves  through  the  wave  equation  and  the  development  of  the 
potential  function  which  is  the  basis  for  the  determination  of  the  fields  in  the  waveguide. 
It  also  contains  an  overview  of  the  field  development  for  the  partially  filled  region  and 
discusses  previous  research  efforts  which  focused  on  the  correction  of  air  gaps  for  the 
waveguide  TR  method.  Chapter  3  highlights  the  development  of  the  Full  Wave  Modal 
Method  for  the  partially  filled  coaxial  line  and  discusses  the  experimental  setup  scenarios 
for  material  parameter  extraction.  This  chapter  also  briefly  discusses  other  methods  for 
solving  the  low  frequency  problem.  Chapter  4  presents  the  experimental  results  and 
discusses  the  potential  sources  of  errors  resulting  from  this  method.  Chapter  5  concludes 
the  thesis  with  a  summary  of  the  effort  and  recommendations  for  future  research. 
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//.  BACKGROUND 


This  chapter  contains  the  background  for  the  characterization  of  thin  samples  in  a 
TEM  coaxial  transmission  line  using  the  Full  Wave  Modal  Method  developed  in  Chapter 
3.  The  terms  waveguide  and  transmission  line  are  used  interchangeably  throughout  the 
document.  Wrapping  a  thin  sample  of  material  around  either  the  inner  or  outer  conductor 
of  a  coaxial  line  presents  a  unique  challenge  to  account  for  the  air-gap  errors  introduced 
by  the  excitation  of  higher  order  modes  in  the  system.  Previous  research  efforts  focused 
on  solving  the  air  gap  problem  using  varying  methods  and  this  thesis  focuses  on  the 
application  of  the  modal  method  to  very  thin  samples  of  high-loss  magnetic  material  in 
the  TEM  coaxial  line  shown  in  Figure  2.1. 

2.1  Transmission/Reflection  (TR)  Method  and  Coaxial  Waveguides 

The  Coaxial  Air- Fine  method  is  one  of  the  most  widely  used  means  to  determine 
the  permittivity  and  permeability  of  a  sample  material.  In  this  method,  a  sample  is 
inserted  between  the  inner  and  outer  conductors.  Coaxial  transmission  lines  come  in  a 
variety  of  sizes  which  are  usually  precisely  designed  to  have  a  characteristic  impedance 
of50Q  .  This  impedance  is  determined  from  the  formula 


where  is  the  relative  dielectric  permittivity  of  the  insulating  material  between  the  inner 

and  outer  conductors  (usually  air)  and  a  and  c  are  the  respective  radii  of  the  inner  and 
outer  conductors.  Most  Network  Analyzers  have  a  characteristic  impedance  of50Q ,  so 
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most  transmission  lines  are  designed  to  mateh  this  value.  In  theory,  a  coaxial 
transmission  line  has  a  lower  operating  frequency  range  of  0  Hz  giving  it  a  wide 
operating  frequency  range  and  making  it  an  ideal  method  for  low  frequency 
measurements.  Typically  the  coaxial  lines  are  named  based  on  their  outer  diameter  and 
Table  2.1  gives  examples  of  several  50Q  coaxial  lines  and  their  approximate  working 
frequency  range  [5].  The  14mm  coaxial  line  used  in  this  research  had  inner  and  outer 
conductor  radii  of  0.122”  and  0.281”  respectively. 

The  TR  method  of  material  characterization  is  widely  employed  because  it 
produces  closed-form  solutions  for  the  determination  of  complex  permittivity  and 
permeability  which  are  based  on  the  measured  scattering  parameters  of  the  sample.  In  a 
coaxial  line,  an  incident  (usually  TEM)  wave  is  introduced  on  the  system  by  a  Network 
Analyzer  which  interacts  with  a  sample  that  fully  fills  the  cross-sectional  area  of  the 
sample  holder  as  shown  in  Figure  2.2  (a).  The  sample  will  reflect  a  portion  of  the  field 
back  to  Port  1  of  the  NWA  which  in  turn  measures  this  reflected  energy  and  calculates 
the  scattering  parameter  S'jj ,  also  known  as  the  reflection  coefficient.  The  remaining 
portion  of  the  energy  will  be  either  absorbed  and/or  transmitted  through  the  sample  and 
the  NWA  measures  this  at  Port  2  and  calculates  5'2i  or  the  transmission  coefficient. 

These  two  values  along  with  the  known  sample  length  are  then  used  in  the  NRW 
algorithm  to  calculate  the  complex  permittivity  and  permeability  of  the  sample  using  the 
formulas  from  [15-16].  Uncertainties  of  course  can  exist  in  this  algorithm,  but  these  have 
been  investigated  by  multiple  researchers  and  can  be  generally  expressed  as 
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where  e-  is  either  the  amplitude  of  S'jj ,  phase  of  S'jj ,  amplitude  of  S 
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length  of  sample  D.  is  the  uncertainty  of  a  where  the  parameter  a 


and  imaginary  parts  of  permittivity  and  permeability  [5]. 


(2.1) 


,  phase  of  ^21  ,  and 
refers  to  the  real 


2.2  The  TEM  Coaxial  Waveguide 

Typically,  a  sample  of  material  is  machined  to  fully  fill  the  cross-section  of  the 
waveguide  so  that  the  NRW  algorithm  may  be  utilized  to  extract  the  material  parameters. 
In  many  cases,  discontinuities  such  as  an  air  gap  like  the  one  shown  in  Figure  2.2  (b)  are 
introduced  to  the  system.  Air  gaps  commonly  present  themselves  as  a  result  of  factors 


Table  2.1:  Examples  of  different  coaxial  transmission  lines  and  their  working  frequency 
range  [5].  The  coaxial  lines  are  typically  referenced  according  to  their  outer  diameter. 
For  this  research,  a  14mm  coaxial  line  was  used. 


Outer  Diameter 
(mm) 

Working  Frequency 
Range  (GHz) 

3.5 

0-34.5 

7.0 

0-18.2 

14.0 

0-8.6 
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Figure  2. 1 ;  14mm  Coaxial  Transmission  Line 


(a) 


t  2a 


Figure  2.2:  A  eoaxial  transmission  line  eontaining  a  sample  that  (a)  fully  fills  the  cross- 
seetional  area  and  (b)  partially  fills  the  eross-seetional  area  with  a  large  outer  airgap. 
Seenario  (b)  is  the  seenario  eonsidered  in  this  researeh. 
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such  as  improper  machining  of  samples,  thermal  contraction  and  expansion,  and  sample 
wear  after  repeated  press  fits  over  time.  Typically  these  gaps  are  small  and  unintentional, 
so  their  effects  are  well  documented.  Air  gaps  are  problematic  because  they  change  the 
fundamental  field  structure  within  the  waveguide  which  can  dramatically  affect 
parameter  extraction  techniques  such  as  the  NRW  algorithm.  To  fully  understand  the 
effects  of  these  air  gaps,  a  thorough  knowledge  of  the  fields  in  the  guiding  structure  is 
required.  Appendix  A  summarizes  the  development  of  the  potential  function,  beginning 
with  Maxwell’s  Equations,  which  is  then  used  to  define  the  fields  in  the  line.  Appendix 
B  shows  the  development  of  the  TEM  and  TM^  fields  in  the  coaxial  waveguide  which  is 
fully  filled  by  free-space.  In  the  case  of  the  empty  or  fully-filled  guide  considered  in 
Appendix  B,  the  TEM  mode  is  the  dominant  propagating  mode  in  the  system.  Higher 
order  TE  and  TM  modes  can  exist  but  are  typically  highly  evanescent  and  considered 
negligible  in  these  cases. 

In  the  partially  filled  guide  as  shown  in  Eigure  2.3,  a  TEM  mode  can  only 
propagate  in  the  partially  filled  region  if  the  z  -directed  phase  front  remains  continuous 
across  the  radial  material  to  air  gap  boundary.  This  implies  that 

Q-jYair-^  _  Q-jYmat^  (2.2) 

where  and  y^^^  are  the  propagation  constants  of  the  wave  through  air  and  the 
material.  In  a  coaxial  line,  the  expression  for  the  z  -directed  propagation  constant  is 

(2.3) 
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where  k^—O  for  the  TEM  mode.  This  implies  that  for  the  TEM  mode,  y  =  k  which  is 

the  frequency  dependent  wavenumber  of  the  material.  In  order  for  (2.2)  to  hold  for  the 
TEM  mode, 

fair  =  Ymat  (2-4) 

which  requires  =  k^^^ .  This  means  that  the  material  must  either  be  free  space  ( 

=  1 )  or  have  relative  permittivity  and  permeability  such  that  =  1  /  //^ .  The 

former  case  is  impossible  with  a  material  other  than  air  and  the  latter  case  is  highly 
unlikely  in  general;  therefore,  a  TEM  mode  physically  will  not  exist  in  the  partially  filled 
region. 


Eigure  2.3:  The  coaxial  transmission  line  containing  a  finite  length  sample  of  material 
which  only  partially  fills  the  radial  cross  section  of  the  guide 
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Since  the  TEM  mode  does  not  exist  in  the  partially  filled  region,  other  higher 
order  TE  or  TM  modes  must  exist  and  propagate  through  the  region.  When  the  NWA 
measures  the  scattering  parameters,  only  the  dominant  mode  is  measured.  Due  to  the  law 
of  conservation  of  energy,  the  energy  that  is  measured  at  the  end  of  the  elosed  system 
must  be  equal  to  that  whieh  is  introduced  to  the  system.  Since  the  TEM  mode  does  not 
propagate  in  Region  II,  the  energy  from  the  propagating  TEM  mode  in  Region  I  must 
eouple  into  the  higher  order  TE  or  TM  modes  at  the  region  boundary.  Since  the  NWA 
measurements  take  into  aceount  only  the  dominant  propagating  mode,  the  energy  coupled 
into  the  higher  order  modes  in  Region  II  is  not  theoretically  accounted  for  leading  to 
errors  in  the  extracted  material  parameters. 

The  first  step  to  accounting  for  these  errors  is  to  determine  what  fields  exist  in  the 
partially  filled  region.  Assuming  axial  symmetry,  the  TE^  fields  in  Region  II  would 
eonsist  of  the  E^,H  components.  The  incident  TEM  fields  from  Region  I  consist  of 

the  E p  and  components.  In  order  for  the  TE  modes  to  exist  in  Region  II,  the  TEM 

incident  fields  must  excite  an  orthogonal  field  (i.e.  e'^^^  must  excite  ),  whieh  for  an 

isotropie  material  considered  in  this  researeh  is  not  possible.  Therefore,  only  TM  higher 
order  modes  ean  exist  in  the  partially  filled  coaxial  region. 

The  forward  propagating  TM^  fields  in  a  general  coaxial  waveguide  are  developed 
in  Appendix  C.  Eor  the  partially  filled  region,  the  z  -direeted  propagation  eonstant 
remains  the  same  in  both  materials,  but  the  radial  components  of  each  region  are  unique. 
This  fact  is  what  allows  for  the  extraetion  of  the  material  properties  of  eaeh  region  from 
the  determination  of  the  propagation  constant.  The  axially  invariant  TM  fields  for  a 
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region  containing  an  outer  air  gap  after  PEC  boundary  conditions  are  applied  at  the  inner 
and  outer  conductors  are 
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where  Yn  -  Jk^  for  the  n*  mode  and 

Vm[Gcp)  =  [jm[Gcp)YQ[aa)-jQ[aa)Y^[ap)] 

Wm  {ocp)  =  [_J„r  {ocp)  Yo  («c)  -  Jo  («c)7„  {ap)\. 

The  next  boundary  condition  which  must  be  enforced  is  the  continuity  of 
tangential  electric  and  magnetic  fields  across  the  material  interface  in  the  region.  When 
evaluated  at  the  material  to  air  interface  at  p  =  h  ,  the  tangential  E  and  H  fields  in  each 
region  are  equal.  This  implies  that 

EAb-)  =  EAP) 

H^{b-)  =  H^{P). 
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Solving  both  equations  for  the  eonstant  B2  o„  yields 
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which  when  equated  produces  the  characteristic  equation  for 

KiAV,  (k,2.b)  W„  -  k^2AV„  (*,,,6)  W,  (k^,„b)  =  0  (2.6) 

where 

kp\n  =  ^fn+kl  and  kp2n  = 

Equation  (2.6)  must  be  solved  for  its  complex  roots,  of  which  there  are  an  infinite 
number.  This  cannot  be  solved  analytically;  therefore,  a  numerical  method  is  needed  to 
find  the  roots.  Methods  such  as  the  Gauss-Newton  or  Levenberg-Marquardt  methods  are 
commonly  used  algorithms  to  solve  equations  of  this  type.  Care  must  be  taken  to  ensure 
that  the  solutions  returned  correspond  to  the  correct  mode,  since  for  magnetic  materials 
the  wavenumbers  are  not  necessarily  sequential  (i.e.  decreasing  complex  components 
with  increasing  real  components). 

2.3  Air  Gap  Correction  Efforts 

Closed  form  solutions  such  as  the  NRW  algorithm  are  the  preferred  method  to 
recover  a  material’s  permittivity  and  permeability  due  to  their  simplicity.  As  discussed 
previously,  these  solutions  have  significant  drawbacks  when  being  used  in  environments 
which  contain  higher  order  modes.  Significant  prior  research  has  been  conducted 
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detailing  attempts  to  derive  air  gap  correetion  factors  for  use  with  methods  such  as  NRW. 
A  few  of  these  methods  include  the  Lumped  Capacitor  and  Inductor  model  developed  by 
the  National  Institute  of  Standards  and  Technology  (NIST)  [I]  and  the  full  wave  modal 
method  performed  in  [8]. 


2.3.1  Circuit  Model  Corrections.  The  correction  factors  developed  by  NIST 
assume  a  single  dominant  mode  in  the  system,  which  for  very  small  air  gaps  is  a 
relatively  good  assumption.  For  the  correction  of  permittivity,  the  partially  filled  line  is 
treated  as  a  system  of  capacitors  in  series  using  the  dimensions  for  the  general  case  of  an 
inner  and  outer  air  gap  shown  in  Figure  2.4  so  that 


J__J_  J_  J_ 


(2.7) 


Figure  2.4:  Partially  Filled  Coaxial  Line  with  Inner  and  Outer  Air  Gaps 
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The  electric  field  distribution  for  a  coaxial  line  is  given  by 


In  —  p 
\a  ) 


where  a  and  c  correspond  to  the  inner  and  outer  conductor  radius.  The  voltage  between 
the  conductors  is  rigorously  defined  as 

(+)  _ 

V  =  -j  E-dl. 

(-) 

and  the  capacitance  of  the  coaxial  line  with  length  L  is  given  by 
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where  R  represents  the  appropriate  radii  from  Figure  2.4.  The  system  of  capacitors 
described  in  (2.7)  is  defined  as 
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where  are  the  corrected  and  measured  real  part  of  the  sample’s  permittivity  and  4 

is  the  real  part  of  the  permittivity  of  air.  This  yields  the  following  corrections  for  the  real 
and  imaginary  parts  of  the  sample’s  permittivity: 
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The  permeability  correction  is  obtained  by  modeling  the  line  as  a  series  of 
inductors  such  that 


L  =L  +L  ■  .  (2.11) 

After  analyzing  the  flux  and  solving  for  the  inductance,  the  corrections  for  the  real  and 
imaginary  parts  of  permeability  are 

_  Ami;  A  ~  A 

McR  ~  j 

^  (2.12) 

Typically,  the  corrected  values  for  permeability  are  much  less  accurate  than  those 
for  the  permittivity  since  H([,  is  continuous  across  the  material  to  air  boundary  while  Ep  is 
not.  The  permittivity  corrections  in  (2.9)  also  break  down  when  in  the  case  of  large  air 
gaps  ( >  (L3/L1)).  While  providing  somewhat  more  accurate  values  in  many  cases, 

these  correction  factors  are  still  very  simplistic  and  based  exclusively  on  the  physical 
dimensions  of  the  line  and  air  gaps.  The  results  of  the  uncorrected  NRW  method  and  the 
NIST  corrected  method  for  the  unknown  dielectric  and  WaveX  materials  are  shown  in 
Figures  2.5  and  2.6  respectively.  As  expected,  the  results  returned  by  the  NRW 
algorithm  for  both  materials  are  poor.  The  corrected  value  for  e\  of  the  dielectric 
material  is  overestimated  but  tends  to  follow  the  general  trend  of  the  fully  filled  scenario. 
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The  corrected  values  of  the  WaveX  measurements  appear  to  correct  the  permeability 
closer  towards  the  actual  values.  The  corrected  permittivity  values  are  quite  poor  until 
after  1.5  GHz.  A  more  exact  solution  requires  consideration  of  the  contributions  from 
these  higher  order  propagating  modes  present  within  the  line.  Some  efforts  approach  the 
problem  using  perturbation  techniques  assuming  a  slight  change  in  the  fields  of  an 
otherwise  fully  filled  line  [10].  These  perturbation  techniques  can  be  good  estimates  but 
still  lack  the  rigor  of  other  more  exact  methods.  Other  methods  use  correction  formulas 
such  as  those  examined  in  [4].  A  rigorous  method  that  has  shown  some  success  for 
dielectrics  and  magnetic  materials  with  air  gaps  is  the  modal  analysis  or  mode-matching 
technique. 

2.3.2  Modal  Analysis.  The  modal  analysis  method  is  a  direct  method  which 
conforms  close  to  physical  reality  by  choosing  the  amplitudes  of  normal  modes  to  satisfy 
boundary  conditions  at  a  discontinuity.  Wexler  [17]  was  one  of  the  early  pioneers  of  this 
method  which  allows  for  the  interaction  effects  of  dominant  and  higher-order  modes 
between  discontinuities  in  a  rectangular  waveguide  to  study  finite  length  obstructions. 
Modal  analysis  considers  an  infinite  number  of  higher  order  modes,  truncated  to  a  finite 
number  for  computational  purposes,  which  are  then  summed  together  and  lead  to  the 
formation  of  the  readily  solvable  linear  algebra  equation 

Ax  =  b. 

The  matrix  A  contains  the  coupling  between  field  modes,  x  is  a  vector  of  unknowns 
which  contains  the  theoretical  reflection  and  transmission  coefficients,  and  b  is  the 
forcing  vector  which  contains  the  known  values  from  the  coupling  of  the  incident  fields 
from  Region  I.  The  reflection  and  transmission  coefficients  correspond  to  the  theoretical 
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S{*J^  and  respectively.  These  values  are  compared  to  the  experimentally  measured 
scattering  parameters  in  the  form  of 


\<5 


(2.13) 


where  S  is  the  specified  tolerance.  A  root  search  algorithm  is  used  to  iterate  the  values 
of  permittivity  and  permeability  in  the  theoretical  calculation  until  (2.13)  is  true  within 
the  tolerance  S . 

A  large  body  of  research  has  focused  on  this  method,  extending  its  applications 
to  many  different  materials  and  environments  [9,  18].  One  such  effort  upon  which  this 
research  is  based  is  the  application  of  the  modal  analysis  technique  to  a  partially  filled 
coaxial  transmission  line. 

In  [8],  Fehlen  investigated  the  scenario  of  a  TEM  coaxial  line  partially  filled  with 
a  linear,  isotropic  and  homogenous  material  using  the  mode-matching  approach.  Both 
non-magnetic  and  magnetic  materials  were  investigated  for  inner  and  outer  air-gap 
scenarios  considering  a  varying  number  of  higher  order  modes.  In  the  freespace  filled 
coaxial  line  used  in  his  work,  he  found  that  32%  of  the  electric  field  component  of  a  TEM 
wave  was  concentrated  in  the  first  2.54mm  extending  radially  from  the  center  conductor 
while  only  17%  of  the  potential  was  contained  in  the  outer  2.54mm.  This  concludes  that 
an  inner  air  gap  prevents  a  large  portion  of  the  potential  from  interrogating  the  sample. 
The  inner  air  gap  likely  leads  to  increased  errors  in  the  measurements,  especially  in  the 
case  of  a  very  large  radial  air  gap.  As  expected,  the  modal  method  produced  more 
accurate  results  when  utilized  with  the  outer  airgap  due  to  the  concentration  of  the 
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electric  field  around  the  center  conductor.  Also,  the  inner  air  gap  is  more  likely  to  induce 
a  drooping  effect  where  the  sample  loses  its  axial  symmetry  due  to  the  lack  of  support 
from  the  center  conductor.  This  elimination  of  symmetry  precludes  the  use  of  the 
corrections  in  (2.9)  and  (2.12)  [5].  The  modal  analysis  technique  returns  good  results  for 
both  materials  although  the  inner  air  gap  scenario  proved  to  be  less  accurate,  particularly 
for  a  magnetic  material.  The  results  for  the  MagRAM  sample  with  a  9  mil  inner  or  outer 
gap  considering  10  modes  is  shown  in  Figure  2.7.  As  expected,  the  outer  gap  scenario 
with  10  higher  order  modes  returned  results  that  agree  well  with  the  fully  filled  case.  The 
detailed  development  of  the  coaxial  line  modal  method  is  summarized  in  Chapter  3. 
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Figure  2.5:  Extracted  complex  and  ju^.  of  an  unknown  dielectric  sample  with  a  123.5 
mil  outer  air  gap  using  the  NRW  algorithm  and  NIST  correction  formulas 
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NRW  - NIST  Correction 


Figure  2.6;  Extracted  complex  and  //,.  of  a  WaveX  sample  with  139  mil  outer  as  gap 
using  the  NRW  algorithm  and  NIST  correction  formulas.  The  sample  is  wrapped  around 
the  inner  conductor  and  radially  fills  approximately  12.5%  of  the  space  between  the  inner 
and  outer  conductors.  The  singular  behaviors  at  1.48  GHz  and  2.95  GHz  is  assumed  to 
be  an  error  in  calibration.  The  corrected  jd^  values  appear  reasonable  across  the  range 

while  the  corrected  values  don’t  appear  reasonable  until  after  1.5  GHz. 
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Figure  2.7:  Measured  complex  permittivity  and  permeability  of  MagRAM  samples  as 
recorded  in  [8].  The  modal  method  with  10  considered  modes  was  used  to  provide 
accurate  corrections  for  both  the  permittivity  and  permeability  of  the  MagRAM. 


2.4  Numerical  Root  Searching 

The  mode-matching  technique  is  a  numerical  technique  which  relies  on  root 
searching  algorithms  to  determine  the  appropriate  values  of  permittivity  and  permeability 
to  satisfy  the  relations  in  (2.13).  There  are  many  different  root  searching  techniques 
which  may  be  used;  however,  two  of  the  more  widely  used  methods  are  the  Gauss- 
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Newton  and  Levenberg-Marquardt  (LM)  methods  for  the  solution  of  non-linear  least 
squares.  The  advantage  to  the  Gauss-Newton  method  is  that  its  rate  of  eonvergence  ean 
approaeh  quadratic;  however,  it  requires  a  relatively  accurate  initial  guess.  A  poor  initial 
guess  can  result  in  failure  for  the  method  to  converge  to  the  actual  solution. 

The  benefit  of  the  LM  method  is  that  it  ensures  the  solution  will  eventually 
converge,  regardless  of  the  accuracy  of  the  initial  guess.  It  converges  because  when  the 
initial  solution  is  far  away  from  the  correct  solution  it  behaves  like  the  method  of  steepest 
descent  which  ensures  a  slow,  but  guaranteed  convergence.  As  the  method  approaches 
the  correct  solution,  it  behaves  like  the  Gauss-Newton  method  which  quickly  converges 
to  the  final  solution.  The  LM  method  is  slower  than  Gauss-Newton,  but  more  robust  and 
is  a  good  choice  when  an  accurate  initial  guess  is  not  known  [12,  14].  In  [8],  well 
quantified  samples  of  Rexolite  and  MagRAM  were  measured,  so  a  combination  of 
Gauss-Newton  and  Muller  complex  root  search  algorithms  were  used.  The  WaveX 
material  investigated  in  this  research  is  not  well  quantified  leading  to  the  selection  of  the 
LM  method  due  to  the  lack  of  an  accurate  initial  guess. 

2.5  Summary 

Observations  from  previous  research  show  that  the  lumped  circuit  model 
corrections  defined  in  [1]  perform  reasonably  well  for  dielectric  materials  but  quite 
poorly  for  heavily  loaded  magnetic  material  samples  in  the  partially  filled  region. 
Additionally,  these  corrections  depend  primarily  on  the  size  of  the  airgap,  becoming 
increasingly  less  effective  as  the  gap  size  increases.  The  mode-matching  technique  used 
in  [8]  returned  accurate  results  for  non-magnetic  materials  filling  only  30%  of  the  cross 
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sectional  area,  especially  when  considering  an  outer  airgap  scenario.  Good  results  for 
high-dielectric  constant  magnetic  material  were  also  obtained  using  this  method.  The 
focus  of  this  research  is  to  extend  this  partially  filled  coaxial  line  modal  analysis 
technique  to  large  radial  air  gap  scenarios  containing  samples  of  the  WaveX  material. 
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III.  METHODOLOGY 


The  partially  filled  TEM  eoaxial  transmission  line  exeites  higher  order  modes  in 
the  regions  eontaining  air  gaps,  requiring  a  method  to  account  for  the  energy  coupled  into 
these  higher  order  modes.  The  full  wave  modal  method  for  this  scenario  was  previously 
developed  and  applied  in  [8]  to  a  system  containing  small  air  gaps  around  either  the  inner 
or  outer  conductor.  This  method  obtained  theoretical  values  for  the  scattering  parameters 

and  S21  for  a  measured  air  gap  size  which  were  then  compared  to  the  measured 


values  5'“'’  and  using 


<s 

<5 


(3.1) 


to  extract  the  permittivity  and  permeability  of  the  material  sample.  This  method  was 
experimentally  verified  to  produce  accurate  results  for  both  a  dielectric  (Rexolite)  and 
magnetic  (MagRAM)  material  after  10  modes.  This  section  contains  a  summary  of  the 
development  of  the  TEM  coaxial  line  modal  method  which  describes  the  relationship 
between  the  fields  within  the  three  regions  of  figure  3.1.  It  also  highlights  the  unique 
factors  considered  by  this  research  for  very  thin  material  samples. 
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Figure  3.1;  The  dominant  TEM  mode  from  Region  1  eouples  into  higher  order  modes  in 
Region  II  excited  by  the  air  gap  discontinuity.  These  Region  2  modes  couple  back  into  a 
TEM  mode  which  propagates  in  Region  3  to  the  Port  2  detector  of  the  NWA. 


3.1  Modal  Method  Development  for  TEM  Coaxial  Transmission  Line 

There  are  three  main  steps  to  develop  the  full  wave  modal  method.  These  three  steps  are; 

1 .  Determine  the  potential  modes  in  the  system  and  define  the  fields, 

2.  apply  boundary  conditions  at  the  interfaces  to  ensure  a  unique  solution, 

3.  and  select  and  apply  a  testing  operator  to  generate  a  linear  system  of  equations. 
As  discussed  in  Chapter  2,  the  partially  filled  TEM  coaxial  transmission  line  contains 
both  a  dominant  TEM  mode  and  higher  order  TM^  modes.  The  TEM  mode  only  exists  in 
the  free-space  filled  Regions  I  and  III  whereas  the  higher  order  TM^  modes  exist  in  all 
three  regions.  The  TEM  and  TM^  field  summaries  for  Regions  I  and  III  are  developed  in 
Appendix  C,  and  considering  only  the  forward  propagating  fields  are  summarized  as 
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where  m  is  either  0  or  1 .  The  charaeteristic  equation  neeessary  to  solve  for  for  the 


TMqj^  modes  is 


Jo  [kpnc) Yq  [kp„a)-J^ [kp„a) Y^  [kp„c) 


=  0 


(3.4) 


where 


kpn  =  yjko  +Yn^  ■ 

The  TM^  fields  in  Region  II  were  developed  in  Chapter  2  and  considering  only 
the  forward  traveling  fields  are 
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and  m  equals  0  or  1 .  The  eharaeteristie  equation  to  solve  for  is 

pnb)  ITo  (t„„6)  -  =  0 

where 


(3.5) 
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(3.7) 


kpin  =  ylYn+ko  and  k^^n  =  ylYn +k2  ■ 

Onee  the  fields  are  identified,  it  is  helpful  to  note  that  only  the  field  components 
which  are  transverse  to  the  direction  of  propagation  will  be  utilized  in  the  modal  analysis 
This  is  due  to  the  application  of  boundary  conditions  at  the  region  boundaries  along  the 
direction  of  propagation.  In  the  case  of  the  z  directed  coaxial  line,  this  means  that  only 
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the  Ep  and  fields  are  used  from  this  point  forward.  The  total  transverse  fields  at  any 
point  in  the  system  are  deseribed  as  the  superposition  of  each  potential  mode  of  the  form 


n=\ 


(3.8) 


n=\ 

where  a*  represents  the  complex  amplitudes  associated  with  each  mode  and  and 
represent  the  respective  transverse  electric  and  magnetic  field  components.  Assuming  the 
material  sample  is  a  simple  media  (i.e.  linear,  isotropic,  and  homogenous)  means  that 

each  value  of  is  constant.  The  propagation  of  the  fields  are  defined  by  the  term 


with  a  propagating  wave  of  the  form  of  in  the  +z  direction. 

Using  the  forms  from  (3.8),  the  modal  equations  for  each  of  the  three  regions 
must  be  developed.  Region  I  contains  a  single  forward  traveling  TEM  mode  which  upon 
scattering  from  the  boundary  at  z=0  will  return  a  reflected  TEM  mode  plus  an  infinite 
number  of  higher  order  TM^  modes.  These  fields  are  represented  as 


z 

17=1 

00 


z  <  0 


(3.9) 


77=1  J 

where  the  n=l  mode  represents  the  dominant  TEM  mode  and  n>l  represents  the  higher 
order  TM  modes. 

As  previously  shown  in  Chapter  2,  a  TEM  mode  cannot  physically  propagate  in 
Region  II;  therefore,  the  incident  TEM  field  from  Region  I  excites  an  infinite  number  of 
forward  traveling  TM^  modes  which  will  then  be  reflected  at  the  z=E  boundary.  This 
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representation  contains  field  vectors  and  amplitudes  along  with  a  unique  propagation 
constant,  all  of  which  are  different  from  those  in  Region  I 


n=\  n=\ 

00  00 


n=\ 


n=\ 


0  <  z  <  Z. 


(3.10) 


Finally,  Region  III  contains  only  forward  propagating  modes  consisting  of  the 
dominant  TEM  mode  plus  an  infinite  number  of  higher  order  TM^  modes  which  are 
represented  as 


'  r+p 

,'^n 


n=\ 


>z>  L. 


(3.11) 


n=\  J 

This  definition  contains  a  phase  shift  of  to  reference  the  fields  to  the  z=L  boundary 
between  Regions  II  and  III  which  simplifies  later  calculations. 

Once  the  fields  have  been  identified,  the  next  step  in  the  process  is  to  apply  the 
boundary  conditions  at  region  interfaces  to  ensure  a  unique  solution  is  maintained.  The 
continuity  of  tangential  fields  in  Region  II  at  the  material  to  free-space  boundary  was 
already  considered  in  the  field  development  summarized  in  (3.5).  Next,  the  continuity  of 
tangential  fields  at  the  z=0  and  z=L  boundaries  must  be  maintained.  This  requires  that 
the  transverse  fields  described  in  (3.8  -  3.10)  must  be  appropriately  equated  and 
evaluated  at  the  appropriate  value  of  z.  This  condition  produces  the  expressions 
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which  is  a  set  of  four  equations  neeessary  to  solve  for  the  unknown  refleetion  and 
transmission  eoeffieients  in  the  system.  These  eoeffieients  found  by  normalizing  the 

amplitudes  of  the  fields  in  eaeh  region  to  the  ineident  TEM  amplitude  and  ean  be 
deseribed  as 
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where  and  eorrespond  respeetively  to  Sjj  and  5*21  of  the  system.  Filling  these  into 
the  system  of  equations  in  (3.12),  rearranging  and  truneating  the  series  to  a  finite  value  N 
yields 
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which  is  now  a  system  of  4  equations  neeessary  with  4N  unknowns.  Considering  only 
the  n=l  mode,  this  system  is  well-posed  and  solvable  for  the  transmission  and  reflection 
coefficients;  however,  the  effect  of  higher  order  (n>l)  modes  cannot  be  considered  by 
(3.14). 

In  order  to  eonsider  the  effects  of  the  higher  order  modes  on  the  transmission  and 
reflection  coefficients,  the  system  of  equations  in  (3.14)  must  be  expanded  to  include  4N 
equations  with  4N  unknowns.  To  aecomplish  this,  a  testing  operator  which  does  not  alter 
the  parameter  definitions  must  be  applied  to  the  system  of  equations.  As  noted  by 
Fehlen,  the  Region  I  fields  serve  as  ideal  testing  operators  since  the  coupling  between  the 
field  modes  is  determined  by  multiplying  the  fields  together  and  integrating  over  the 
eross-section  of  the  coaxial  line.  The  scattering  parameters  remain  constant  over  the 
surface  integral,  so  the  testing  operator  does  not  impaet  these  parameters.  A  benefit  of 
choosing  the  Region  I  fields  for  the  testing  operator  is  the  orthogonality  of  fields  in  a 
fully  free-space  filled  region  sueh  as  Regions  I  and  III.  This  means  that  ^ox  m^n  the 
fields  do  not  couple  into  eaeh  other,  reducing  these  integrals  into  a  set  of  simplified 
diagonal  matriees.  The  Region  I  and  II  fields  are  not  orthogonal;  therefore,  a  full  set  of 
4N  X  4N  matriees  will  be  formed  accounting  for  all  higher  order  mode  eoupling.  The 
testing  operator  is  applied  to  each  term  in  (3.14)  sueh  as 


cs 

\K-{]  (m  =  l,2,...,A), 


(3.15) 


cs 

choosing  the  appropriate  operator  for  the  electric  or  magnetic  field.  Applying  the  testing 
operator  to  (3.14)  yields 
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(3.16) 
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which  describe  the  eoupling  between  the  field  eoupling  between  regions.  A  separate 
integration  to  deseribe  the  magnetie  field  eoupling  between  the  Region  1  and  111  fields  is 
not  neeessary  due  to  the  relation  of  the  eleetrie  and  magnetie  fields 


where  is  the  wave  impedance  for  the  nth  mode.  This  relationship  simplifies  the 
produet  and  integration  of  the  testing  operator  fields  to 
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whieh  is  simply 
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Using  the  testing  operator  definitions  from  (3.16)  and  (3.17),  the  system  of  equations  in 
(3.14)  is  rewritten  as 
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where  eaeh  element  in  the  matrix  is  an  NxN  sized  submatrix.  Equation  (3.18)  is  now  a 
well  posed  linear  system  of  equations  readily  solvable  by  linear  algebra  teehniques. 

The  next  step  is  to  summarize  the  transverse  fields  in  eaeh  region  whieh  are  used 
in  the  integral  caleulations  of  the  testing  operators.  The  transverse  fields  for  a  z  direeted 
eoaxial  line  are  the  radial  and  axial  eomponents  of  the  eleetrie  and  magnetie  fields.  In  the 
ease  of  the  particular  fields  for  this  research  summarized  in  (3.2)  and  (3.5),  the  radial 
electric  field  and  axial  magnetic  field  ( )  are  the  transverse  fields.  Normalizing 
these  fields  to  the  incident  TEM  wave  by  the  constant 
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for  free-space  filled  Regions  I  and  III  and 
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for  the  partially  filled  Region  II  where 

K  [ap)  =  [  [ap)  7o  («a)  -  Jo  (aa)  («p)] 

K  M  =  [Jm  («P)  Yx,  («c)  -  Jo  («c)  («p)] , 

and  are  determined  by  the  characteristic  equations  (3.4)  and  (3.7),  and  subregions  I 
and  2  respectively  correspond  to  the  air-filled  and  material-filled  portions  of  Region  2. 

Using  the  fields  described  above,  the  testing  operator  integral  calculations  may  be 
performed.  The  integral  across  the  cross  section  of  a  coaxial  line  involves  integrating 
along  both  the  axial  and  radial  components.  Since  this  research  assumes  axial-invariance, 
the  axial  integral  will  produce  a  factor  of  27i  for  each  case,  which  is  simply  normalized 
out  of  all  the  calculations. 

The  coupling  of  the  electric  fields  in  the  free-space  filled  Regions  I  and  III  is 
defined  from  (3.16)  as 
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Knowing  that  fields  in  these  regions  are  orthogonal  to  eaeh  other,  only  the  self  eoupling 
fields  (m=n)  will  produee  a  value  other  than  zero.  For  the  TEM  (m=n=l)  ease, 
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...Wronskian  relationship  [3]. 

The  mode  eoupling  of  fields  from  Region  I  to  Region  II  are  defined  from  (3.16)  as 
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These  two  integrals  may  be  split  based  upon  the  radial  discontinuity  to  form  two  separate 
integrations  for  each  case  corresponding  to  the  fields  in  the  material  filled  region  and 
those  in  the  air-gap  region.  These  integrals  are  then  summed  to  produce 
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Each  case  must  be  analyzed  for  the  TEM  (m=l)  and  TM  (m>l)  testing  operators  which 
produce  distinctly  different  results.  The  development  of  the  integrals  in  (3.24)  is 
discussed  in  depth  in  [8],  and  result  in  the  following  solutions: 
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(3.28) 


Using  the  results  from  (3.22)  through  (3.28),  the  matrix  deseribed  in  (3.18)  ean  be  filled 
in  and  solved  for  the  unknown  reflection  and  transmission  coefficients  which  correspond 
to  the  theoretical  S'jj  and  5*21  parameters.  Now,  a  root  search  algorithm  may  be  used  to 


solve  (3.1) 

to  return  the  relative  complex  permittivity  and  permeability  of  the  material  under  test.  As 
previously  noted  in  Chapter  1 ,  a  strong  reflection  response  is  not  expected  due  to  the 
small  thickness  of  the  material;  therefore,  the  Sn  values  will  likely  not  provide  a 
sufficient  level  of  fidelity.  To  overcome  this,  two  different  length  material  samples 
should  be  measured  and  the  theoretical  scattering  parameters  for  each  developed  using 
the  modal  method.  Then,  the  S'jj  values  from  (3.1)  could  be  replaced  by  a  new  ^21  value 


independently  derived  from  the  new  length  sample.  The  comparison  equation  now 
becomes 
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which  contains  two  independently  derived  values  of  5*21 . 
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3.2  Other  Low  Frequency  Methods 

In  addition  to  the  coaxial  modal  method,  a  stripline  approach  is  briefly  explored. 
This  method  is  an  attempt  to  develop  an  approach  which  can  take  advantage  of  an 
analytical  technique  such  as  the  NRW  algorithm  to  rapidly  and  accurately  determine  the 
complex  permittivity  and  permeability  of  a  material.  The  stripline  easily  allows  for  thin 
samples  of  material  to  be  investigated  making  it  a  viable  alternative  to  the  coaxial  line. 

One  drawback  of  striplines  is  the  relative  difficulty  in  calibration.  To  overcome 
this,  a  post  collection  calibration  is  accomplished.  Considering  the  geometry  shown  in 
Figure  3.2  where  Regions  A  and  B  correspond  to  the  Port  1  and  Port  2  sides  of  the 
material  sample  of  thickness  d,  the  measured  scattering  parameters  relate  to  the  scattering 
parameters  of  the  material  under  test  by 


^MUT 

*^11  “  ^21^11  ^12 

om  ^MUT  oS 

^22  “  ^12^22  ^21 

om  oMC/r  oS 

^21  “  ^21^21  ^21 

^^12  ”  ^^12^^12  ^12* 


(3.30) 


Rearranging  the  terms  in  (3.30)  gives  the  S-parameters  for  the  material  under  test 


cMUT  _  ‘^11 

‘^21*^12 


I^MUr  _  ^22 

22 

^12^21 


^MUT  _  ^^21 

21  ^  A  ^B 

'^2l'^2l 


cMUT  _  ‘^12 

^^12  -  ^.B  ^.A 

‘^12*^12 


(3.31) 

(3.32) 

(3.33) 

(3.34) 
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where  5"”  is  the  uncalibrated  measured  value  returned  by  the  Network  Analyzer.  The 
denominators  of  (3.31)  -  (3.34)  are  based  upon  measured  results  from  either  an  air  filled 
line  or  a  shorted  line  in  order  to  provide  a  calibration  factor.  When  considering  the  air 
filled  line,  the  sample  region  thickness  d  is  assumed  to  be  the  same  thickness  as  the 
MUT. 

To  solve  for  the  denominators  associated  with  the  5'2i  and  5j2  values,  an  empty 
line  is  analyzed.  Considering  an  arbitrary  “sample-region”  of  thickness  d,  the  measured 
5*21  and  Si2  values  of  air  are 


^21 


’12 


which  when  rearranged  yield 


_  ^m,air  jkd 

^21^21  ”  ^21  ^ 


_  ^m,air  jkd 

^12^12  ”  ^12  ^ 


(3.35) 


The  results  from  (3.35)  are  fdled  back  in  to  (3.33)  and  (3.34)  respectively.  To  solve  for 
the  unknown  denominators  from  (3.31)  and  (3.32)  ,  measurements  using  a  short  are 
required.  First,  a  short  is  inserted  at  d  =  0  and  S-parameters  are  recorded.  The  measured 
of  the  short  is 


since  the  reflection  coefficient  of  a  short  is  simply  -1.  Rearranging  results  in 

52^12  (3-36) 


which  is  used  by  (3.31). 
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The  final  unknown  denominator  ( 5'j2*S'2i )  is  slightly  more  difficult  to  solve  for  due 
to  phase  progression  issues.  Considering  the  short  has  a  thickness  of  L, 

(3.37) 

where  B  is  the  total  distance  from  the  port  1  face  of  the  short  to  port  2.  Knowing  that 

_  o(5— i)  o/, 

^12  “  ^12  ^12 

(3.37)  essentially  amounts  to  a  phase  delay  based  on  the  thickness  L  of  the  short.  This 
means  that 


c(B-i)  _  qB  jkL 


An  identical  process  is  used  to  solve  for  which  yields  an  identical  phase  delay  of 

.  Inserting  these  values  into  (3.37)  and  solving  for  the  resulting  unknown 
yields 


o  m  ,sh 

cB  cB  _  “^22 
^j2kL 


(3.38) 


which  is  the  unknown  for  (3.32). 

Equations  (3.31)  -  (3.34)  are  then  solved  for  the  calibrated  scattering  parameters 
of  the  material  under  test.  These  calibrated  S-parameters  are  then  simply  evaluated  using 
the  NRW  algorithm  for  an  analytical  closed  form  solution  of  the  complex  permittivity 
and  permeability  of  the  material. 
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A - d 


B 


Figure  3.2:  Stripline  eontaining  material  sample  region  of  arbitrary  thiekness  d.  In 
addition  to  the  material  sample,  measurements  are  also  taken  using  both  air  and  a  short  as 
the  “sample”  to  provide  the  neeessary  measurements  to  ealibrate  the  measured  data. 

3.3  Summary 

The  modal  method  for  a  TEM  eoaxial  transmission  line  is  deseribed  in  detail  in 
the  preceding  sections.  The  concept  of  using  multiple  independent  82^  measurements  to 

solve  for  and is  introduced  to  take  the  place  of  using  both  and  5*21  values. 
This  was  developed  based  on  the  thought  that  thin  samples  containing  very  large  air  gaps 
would  have  an  almost  exclusive  transmission  component  with  very  little  scattering  back 
towards  the  incident  port.  A  method  for  using  a  low  frequency  transmission  line  to 
obtain  calibrated  scattering  parameters  for  thin  material  samples  is  also  briefly  discussed. 
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IV.  RESULTS 


The  modal  method  developed  in  the  previous  ehapter  was  used  to  analyze  samples 
whieh  partially  filled  a  TEM  eoaxial  transmission  line  at  low  frequeneies.  This  method 

ealculates  the  theoretieal  scattering  parameters  and  Sfi  which  are  then  compared  to 
the  measured  S-parameters  in  a  ID  or  2D  complex  root  search.  Different  length  samples 
are  measured  to  obtain  two  independent  transmission  measurements.  The  additional  S21 
value  may  be  used  in  place  of  the  ^  parameter  in  the  event  that  the  measured  S'j  j  value 

is  too  small  to  accurately  use  to  perform  the  minimization  in  equation  (3.1).  This  is  a 
distinct  possibility  given  the  large  air  gaps  considered  in  this  research  preclude  much  if 
any  scattering  from  the  material  sample.  Like  in  (3.1),  the  theoretical  values  of 
permittivity  (^)  and  permeability  (//)  are  varied  for  each  frequency  step  until  the 
difference  between  the  theoretical  and  measured  S-parameters  satisfy 


<5 


<S. 


(4.1) 


4.1  Test  Setup 

The  test  setup  consisted  of  a  14mm  coaxial  waveguide,  an  Agilent  E8362B 
Network  Analyzer  and  associated  cables,  and  samples  of  an  unknown  dielectric  material 
and  the  magnetic  material  WaveX  (WX-AOXO).  The  dielectric  samples  are  used  to 
determine  the  effectiveness  of  the  multiple  5*2  j  technique  which  then  is  applied  to  the 
WaveX  samples. 
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The  coaxial  test  fixture  is  shown  broken  apart  in  Figure  4.1  and  fully  assembled  in 
Figure  4.2.  The  waveguide  is  mounted  to  a  plate  with  guiding  rods  to  ensure  proper 
alignment  each  time  the  fixture  was  opened.  The  sample  holder  is  inserted  between  the 
two  ports  and  aligned  via  aligning  pins  and  then  secured  with  three  tension  clamps  to 
ensure  no  energy  leakage  at  the  interface. 

The  inner  conductor  of  the  coaxial  line  has  a  radius  of  122  mils  and  the  outer 
conductor  has  a  radius  of  281  mils.  The  waveguide  has  three  different  length  sample 
holders  (880  mil,  4001  mil,  and  4381  mils)  with  the  first  two  used  in  this  research.  The 
waveguide  was  calibrated  using  a  two  port  TRL  calibration  as  described  in  [7],  using  the 
appropriate  sample  holder  as  the  line  standard.  This  sets  the  measurement  calibration 
planes  at  the  port  1  (z=0)  and  port  2  (z=h)  faces  of  the  sample  holder.  The  sample 
lengths  used  for  both  the  dielectric  and  WaveX  ranged  from  500-1500  mils  which 
sometimes  did  not  completely  fill  the  sample  holder  lengthwise.  The  offsets  from  the 
calibration  planes  were  noted  and  used  to  appropriately  advance  or  delay  the  measured 
scattering  parameters  to  account  for  the  unfilled  length  of  the  sample  holder. 

The  center  conductor  of  the  line  is  held  in  place  by  two  “bullets”  which  are  press 
fit  into  each  end  of  the  center  conductor  and  mate  with  the  corresponding  Portl  and  Port2 
ends  of  the  fixture.  It  was  assumed  that  the  center  conductor  remained  rigidly  fixed  to 
ensure  that  only  axially  symmetric  modes  were  excited  in  the  sample  region.  A 
frequency  sweep  from  50  MHz  to  3.0  GHz  is  conducted  by  collecting  data  at  401  discrete 
points  along  the  band. 

Two  different  sample  lengths  of  dielectric  were  machined  to  a  wall  thickness  of 
35.5  mils,  yielding  a  123.5  mil  outer  air  gap.  These  samples  fdled  approximately  22.3% 
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of  the  radial  surface  area  between  inner  and  outer  conductor.  An  additional  dielectric 
sample  was  machined  to  completely  fill  the  sample  holder  to  determine  the  baseline 
material  parameters  using  the  NRW  algorithm.  This  is  the  value  that  is  used  in 
comparison  to  the  partially  filled  scenarios. 


Figure  4. 1 ;  Coaxial  Test  Fixture  with  Sample  Holder  and  Inner  Conductor  Exposed 


Figure  4.2:  Assembled  Coaxial  Test  Fixture 


The  WaveX  samples  were  cut  from  sheets  of  10  mil  and  20  mil  thickness  and  then 
wrapped  around  the  center  conductor  resulting  in  outer  air  gaps  of  149  mils  and  139  mils 
respectively.  These  samples  respectively  filled  approximately  6.3%  and  12.6%  of  the 
radial  surface  area  between  inner  and  outer  conductor.  The  WaveX  sheets  contain  a  self- 
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adhesive  backing  which  allowed  them  to  be  secured  to  the  center  conductor  without  the 
use  of  any  tapes  or  pastes.  The  samples,  once  secured  to  the  center  conductor,  contained 
slight  imperfections  as  the  samples  were  prepared  and  wrapped  by  hand  without  the  use 
of  precision  machining.  The  most  prevalent  source  of  error  in  this  method  is  likely  to  be 
in  the  axial  plane  as  the  edges  of  the  material  were  not  always  completely  mated.  This  is 
expected  to  have  a  greater  impact  on  the  magnetic  field  as  it  is  (j)  directed  and  will 
experience  a  slight  discontinuity  at  the  boundary.  Additionally,  a  sample  which  fully  fills 
the  sample  holder  was  crafted  by  wrapping  multiple  layers  on  top  of  each  other  since  a 
machined  sample  was  unavailable  for  this  research.  This  yielded  a  crude  sample  for 
which  the  NRW  algorithm  could  be  used  to  extract  the  permittivity  and  permeability  to 
get  a  rough  baseline  for  comparison. 

As  previously  mentioned,  strong  reflection  responses  were  not  anticipated  since 
both  materials  fill  less  than  25%  of  the  radial  surface  area.  Different  length  samples  were 
prepared  to  obtain  multiple  independent  transmission  ( 5*21 )  measurements.  For  the 
dielectric  material,  only  a  ID  complex  root  search  is  required  to  find  since  //^  =  1  is 
known;  therefore,  only  one  5*21  measurement  is  needed.  A  2D  complex  root  search  also 
performed  for  the  dielectric  using  the  measured  and  ^21  values  which  compares 
fairly  well  to  the  fully-filled  and  ID  solution.  Since  the  WaveX  contains  both 
permittivity  and  permeability,  a  2D  complex  root  search  is  needed.  The  magnitudes  of 
the  S-parameters  for  the  WaveX  sample  are  shown  in  Figure  4.4.  Since  the  |  5'j  J  for 

WaveX  is  relatively  large,  a  root  search  using  both  and  S21  values  could  be  used. 
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This  thesis  presents  data  obtained  using  and  5'2i  for  the  root  search  as  well  as  data 


obtained  using  two  independently  obtained  ^21  measurements. 

4.2  Modal  Method  Performance 

The  modal  method  was  used  to  extract  the  relative  permittivity  and  permeability 
from  samples  of  a  dielectric  and  WaveX,  both  of  which  only  partially  fill  the  waveguide 
resulting  in  large  outer  air  gaps.  The  performance  of  the  modal  method  for  both  materials 
is  described  here. 

4.2.1  Performance  for  a  thin  dielectric.  The  modal  method  performed  quite 
well  for  the  unknown  dielectric  samples  measured  in  this  research.  Performing  a  ID  root 
search  using  the  measured  *^21  values,  the  modal  method  considering  20  modes  was  used 
across  the  entire  frequency  range  to  solve  for  the  relative  permittivity.  Figure  4.3  shows 
the  extracted  permittivity  values  from  the  modal  method  compared  to  the  fully  filled 
permittivity  obtained  via  the  NRW  algorithm.  The  modal  method  results  clearly 
converge  to  the  fully  filled  results  after  only  a  few  frequency  steps  and  agree  well  with 
the  fully  fdled  case  across  the  frequency  range  with  a  couple  of  exceptions.  The  fully 
filled  scenario  consistently  exhibits  spikes  at  approximately  1.48  GHz  and  2.95  GHz 
which  are  likely  the  frequencies  that  correspond  to  a  half  and  full  wavelength  for  the 
empty  line  standard  during  calibration. 

A  2D  complex  root  search,  iterating  both^^  and  ,  was  then  conducted  using 
both  the  measured  S'jj  and  5*21  values  from  the  same  dielectric  sample  used  in  the  ID 
case.  The  magnitudes  of  the  measured  S'jj  and  5*21  values  are  shown  in  figure  4.4.  As 
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expected,  the  transmitted  component  greatly  dominates  the  reflected  component  across 
the  entire  frequency  range.  The  extracted  values  of  and  returned  from  the  2D 
search  are  shown  in  Figure  4.5.  In  both  the  ID  and  2D  searches,  a  higher  than  expected 
value  for  s''  was  reported.  This  would  indicate  that  the  material  is  of  a  slightly  more 

lossy  nature  than  originally  anticipated. 

Another  flnding  discovered  for  the  dielectric  samples  is  the  number  of  modes 
required  to  achieve  convergence.  Figure  4.6  clearly  shows  that  the  results  appear  to 
converge  after  the  use  of  only  10  modes.  This  finding  coincides  very  well  with  the 
results  reported  in  [8]  for  samples  of  Rexolite®.  A  perturbation  analysis  was  also 
conducted  using  the  method  described  in  [10]  to  obtain  an  approximate  solution  for  s^. 
This  method  considers  that  the  permittivity  of  the  material  is  only  slightly  perturbed  from 
that  of  free  space  yielding  an  expression  for  Sf^jjj.  of 

^MUT  *  ^0  (1  ) 

721n(^,7"/^,t)ln(c/a) 

{Rq!  Sq)1  ln(£i/a) 

{Rq  I  S(^)l  ln(h/a) 

The  comparison  of  the  perturbation  method  to  the  fully  filled  case  is  shown  in  Figure  4.7. 
It  is  evident  from  this  that  the  higher  order  modes  do  carry  a  significant  enough  amount 
of  energy  that  they  need  to  be  included  in  the  analysis. 
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20  Modes  123.5  mil  Gap - Fully  Filled 


Figure  4.3:  Calculated  of  a  dielectric  with  a  123.5  mil  outer  air  gap  using  20  modes. 

The  spikes  in  the  Fully  Filled  NRW  values  are  the  result  of  an  unidentified  error  in  the 
calibration  kit;  otherwise,  the  calculated  values  agree  very  well  with  the  fully  filled 
values.  These  results  were  obtained  from  a  ID  root  search  using  S21  values. 


Frequency(GHz) 


Frequency(GHz) 


Figure  4.4:  Magnitude  of  and  5*21  for  a  0.5”  long  sample  of  dielectric  with  a  123.5 
mil  outer  air  gap.  As  expected,  |  5'j  J  is  extremely  small  due  to  the  large  air  gap.  These 
magnitudes  are  used  to  perform  a  2D  root  search  to  solve  for  and  for  the  dielectric. 
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20  modes  123.5  mil  Outer  Gap - Fully  Filled 
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Figure  4.5:  Extracted  complex  permittivity  and  permeability  of  a  0.5”  dielectric  sample 
with  a  123.5  mil  outer  air  gap  considering  20  modes  compared  with  the  extracted  values 
from  the  fully  filled  case.  Much  like  the  results  returned  from  the  ID  root  search,  the 
complex  permittivity  agrees  very  well  with  the  fully  filled  results  aside  from  the 
singularities  in  the  fully  filled  case.  As  expected,  the  extracted  permeability  values  are 
//'  « 1  and  //"  «  0. 
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123.5  mil  Outer  Gap - Fully  Filled 


Figure  4.6:  Measured  eomplex  permittivity  at  1  GHz  shown  as  s^-e[-  js".  The 
results  appear  to  converge  after  approximately  10  modes  for  the  low  dielectric  constant 
material  measured. 


Perturbation 


Fully  Filled 


Figure  4.7:  A  perturbation  analysis  is  performed  as  described  in  [10]  to  find  the  complex 
permittivity  s[  of  the  dielectric.  This  assumes  that  the  fields  in  the  sample  region  behave 
like  the  empty  guide  dominant  mode,  only  slightly  perturbed  by  Af  . 
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4.2.2  Performance  for  WaveX.  The  modal  method  performed  quite  poorly 
for  the  samples  of  20  mil  thick  WaveX  material  with  a  139  mil  outer  air  gap.  A 
maximum  of  5  modes  was  used  in  the  modal  analysis  as  the  inclusion  of  additional  higher 
order  modes  caused  the  algorithm  to  break  down.  As  expected,  the  first  mode  excited  in 
the  sample  region  (  IMqq  )  carries  the  vast  majority  of  the  energy  while  the  other  higher 
order  modes  all  evanesce  very  rapidly  as  shown  in  Figure  4.8 

The  original  assumption  was  that  the  magnitude  of  the  returned  S'jj  measurements 
would  be  insufficient  to  use  in  a  root  search,  leading  to  a  need  to  obtain  two  independent 
5*21  measurements.  As  shown  in  Figure  4.8,  the  magnitudes  were  large  enough  to  be 

used,  so  the  parameter  extraction  was  conducted  using  both  methods.  The  results  from 
both  methods  are  shown  in  Figures  4.9  and  4.10. 

The  extracted  /u^  values  surprisingly  yielded  relatively  accurate  results,  but  the 

returned  values  for  were  very  inaccurate,  particularly  at  the  lower  frequency  ranges. 
The  values  agreed  relatively  well  with  the  GTRI  in-plane  permeability  shown  in 
Figure  D.3  for  the  “Old”  sample.  They  also  agree  relatively  well  with  the  fully  filled 
solution  obtained  from  the  crafted  sample  which  fully  filled  the  guide.  The  values 
extracted  had  very  high,  negative  imaginary  parts  and  very  large  real  components,  both  of 
which  are  not  readily  explained.  They  should  have  agreed  with  the  out-of-plane 
permeability  results  shown  in  Figure  D.  1 .  This  phenomena  was  consistent  with  both 
approaches  to  the  root  search. 
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Figure  4.8;  The  dB  magnitude  of  the  measured  Sjj  for  a  1.48”  long  sample  of  WaveX 
with  a  139  mil  outer  air  gap. 


- Fully  Filled  -  5  Modes  139  mil  Outer  Air  Gap 


Figure  4.9;  Extraeted  and  extraeted  from  a  WaveX  sample  with  a  139  mil  outer  air 
gap  eonsidering  4  modes  using  and  S21 .  The  permeability  agrees  reasonably  with  the 
fully  filled  solution  but  the  permittivity  eontains  a  great  deal  of  error. 
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4  Mode  139  mil  Outer  Air  Gap - Fully  Filled 


Figure  4.10;  Extracted  and  //^extracted  from  a  WaveX  sample  with  a  139  mil  outer 
air  gap  considering  4  modes.  The  permeability  agrees  relatively  well  with  the  fully  filled 
solution  but  the  permittivity  contains  a  great  deal  of  error. 

One  finding  discovered  during  the  WaveX  testing  is  the  anisotropy  of  the 
material.  Prior  to  receiving  the  GTRl  results,  a  sample  was  measured  in  a  low  frequency 
stripline  to  obtain  a  baseline  for  comparison  to  the  modal  method.  The  measurements 
were  taken  as  noted  in  Section  3.2  and  the  NRW  algorithm  was  used  to  extract  the  values 
of  and  of  the  material.  Due  to  a  high  level  of  noise  in  the  measurements,  a 

windowing  function  was  used  to  isolate  the  primary  response.  The  data  was  transferred 
to  the  time  domain  on  the  network  analyzer  where  the  primary  response  was  located.  The 
window  was  set  at  the  nulls  of  the  primary  response  and  data  transferred  back  into  the 
frequency  domain.  This  resulted  in  much  cleaner  S-parameter  measurements,  eliminating 
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a  large  portion  of  the  noise.  As  shown  in  Figure  4.11,  the  permeability  returned  a  similar 
trend  to  the  modal  method  results  and  GTRI  results  although  shifted  lower  in  magnitude. 
The  values  of  settled  down  to  a  reasonable  value  after  approximately  500  MHz  and 

then  stayed  within  a  reasonable  range  up  to  3  GHz.  The  values  of  s”  oseillated  widely 
aeross  the  frequency  range  which  is  not  readily  explained. 

As  an  anisotropic  material,  WaveX  will  have  different  values  of  permittivity  and 
permeability  depending  on  whether  the  sample  is  “in-plane”  or  “out-of-plane”  with  the 
respective  electric  and  magnetic  fields.  For  example,  when  measuring  a  sample  wrapped 
around  the  center  conductor  of  the  coaxial  line,  the  “out-of-plane”  permittivity  and  “in¬ 
plane”  permeability  will  be  extracted.  This  is  because  the  radially  directed  electric  field 
travels  out  of  the  plane  of  the  material  while  the  azimuthally  directed  magnetic  field  is  in 
plane  with  the  material.  In  the  stripline,  the  material  is  oriented  differently  so  that  the  z 
directed  electric  field  is  in  plane  with  the  sample  and  the  circulating  magnetic  field  is  also 
in  plane  with  the  sample.  This  explains  the  similarities  and  differences  between  the 
extracted  permittivity  and  permeability  values  from  the  stripline  and  the  coax. 

A  sample  of  plexiglass  was  also  measured  in  the  stripline  using  the  same  method 
to  ensure  that  accurate  results  were  obtained.  The  extracted  permittivity  and  permeability 
values  are  shown  in  Figure  4.12.  Using  a  baseline  value  of  —  2.6  and  /d^—  \  the 
extracted  values  from  the  stripline  yielded  accurate  results  after  approximately  500  MHz. 
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Figure  4.11:  The  complex  permittivity  and  permeability  of  WaveX  calculated  from 
stripline  measurements  are  shown  as  =  a^  -  j^r  l^r  =  l^r~  Jf^r- 


Stripline  - Actual 


Figure  4.12:  Permittivity  and  permeability  of  plexiglass  calculated  from  stripline 
measurements  are  shown  as  a^  =  a^  -  ja^  and  ju^  =  -7//,. . 
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4,3  Wavenumber  Calculation. 


The  accurate  calculation  of  wavenumbers  is  critical  to  the  success  of  the  modal 
method.  The  wavenumbers  are  calculated  by  solving  for  the  roots  of  (3.7)  which  is 
rewritten  here  as 

(*,2,*)  (y„6)  -  ir,(l:^„b)  =  0.  (4.2) 

There  are  an  infinite  number  of  roots  to  this  equation  as  a  result  of  the  oscillatory  nature 
of  the  Bessel  functions  which  must  be  found  in  sequential  order. 

Given  the  relatively  large  outer  air  gaps  considered,  the  initial  approach  to  solving 
for  the  roots  of  (4.2)  was  to  use  a  perturbation  approach,  also  known  mathematically  as  a 
Homotopy  method.  The  method  breaks  the  material  thickness  into  a  number  of  slices 
which  are  slowly  added  to  each  other  until  the  desired  material  height  is  reached.  The 
initial-root  guesses  used  for  the  first  iteration  are  the  roots  to  the  characteristic  equation 
of  the  unfilled  guide.  The  empty  waveguide  is  then  “stepped”  out  by  one  thin  layer  of 
material  and  (4.2)  is  solved  using  the  empty  guide  solution  as  the  initial  guess.  This 
process  is  repeated  using  the  previous  step’s  solution  as  the  initial  guess  for  the  next 
iteration.  This  is  done  until  the  final  thickness  h  of  the  material  is  reached.  The  final 
solutions  at  that  point  are  the  actual  values  of  in  the  sample  region.  This  assumes  that 

the  roots  are  continuously  dependent  on  the  sample  thickness. 

This  method  is  simple  and  effective,  giving  a  high  likelihood  that  the  numerical 
technique  used  to  solve  (4.2)  will  converge  to  the  proper  root.  The  major  drawback  to 
this  lies  in  the  step  size  (S)  required  to  ensure  that  the  solution  converges  to  the  proper 
root.  This  step  size  is  dependent  on  both  the  ultimate  thickness  h  and  the  material 
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properties  of  the  material  under  test.  A  step  size  of  5  =  hl  (15  A),  where  N  is  the  total 


number  of  modes,  was  found  to  suffieiently  ensure  f  eonvergenee  to  the  eorreet  values 
for  materials  with  low  to  medium  values  of  and  [11].  This  value  of  5  works  well 

for  the  dieleetric  material  but  must  be  inereased  for  higher  order  modes  of  the  WaveX. 
Since  the  total  number  of  iterations  is  linear  with  N,  the  runtime  required  by  this 
algorithm  can  quickly  become  very  large.  While  accurate,  this  method  quickly  becomes 
impractical  due  to  its  runtime;  therefore,  a  faster  method  was  needed. 

In  [8],  a  Muller  root  search  method  was  used  since  the  first  roots  for  a  highly 
loaded  material  are  often  in  close  proximity  to  each  other.  Muller’s  method  is  a  variation 
of  the  secant  method,  requiring  three  initial  guesses  to  the  equation.  It  fits  a  parabola 
through  these  three  points  and  examines  the  intersection  of  the  parabola  with  the  complex 
plane.  This  intersection  point  becomes  the  next  initial  guess  and  the  process  is  repeated 
until  the  root  of  the  function  is  located.  Muller’s  method  works  well  with  complex  roots 
provided  the  function  is  analytic.  Since  Bessel  functions  of  the  first  kind  are  analytic 
everywhere  and  Bessel  functions  of  the  second  kind  are  singular  only  at  the  origin, 
Muller’s  method  is  ideally  suited  for  finding  the  roots  of  (4.2).  The  initial  guesses  for 
low  frequency  measurements  of  low-loss  dielectric  materials  are  well  estimated  by  the 
roots  of  the  characteristic  equation  for  the  unfilled  coaxial  line.  Muller’s  method  quickly 
and  accurately  locates  the  roots  of  the  partially  filled  case  which  are  indeed  very  close  to 
those  of  the  empty  waveguide. 

The  wavenumber  calculations  for  heavily  loaded  magnetic  materials  such  as 
WaveX  are  much  more  difficult  to  solve  since  the  solution  is  extremely  sensitive  to  the 
initial  guess.  In  [8],  a  method  using  asymptotic  expansions  of  the  Bessel  functions  was 
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used  to  determine  the  initial  guesses  for  use  with  the  Muller  root  seareh.  Using  this 
approach,  the  correct  wavenumbers  were  quickly  found  for  a  MagRAM  filled  coaxial  line 
containing  a  9  mil  inner  or  outer  air  gap. 

Both  the  Muller  method  and  the  step/perturbation  method  were  implemented  to 
find  the  roots  of  (4.2)  for  the  dielectric  material.  The  differences  between  the  two  results 

were  found  to  be  on  the  order  of  1x10”^  and  the  runtime  of  the  Muller  method  was  much 
faster  than  the  Homotopy  method.  The  comparison  between  the  wavenumbers  for  the 
empty  waveguide  and  the  guide  partially  filled  with  dielectric  is  shown  in  Table  4.1. 

Both  methods  were  evaluated  to  solve  for  the  roots  of  (4.2)  for  the  WaveX  samples.  The 
first  mode  returned  by  both  methods  agreed  well,  but  the  perturbation  method  returned 
several  additional  higher  order  modes  not  found  by  the  Muller  method  search.  The 
comparison  between  the  two  is  shown  in  Table  4.2. 

With  the  disparity  between  the  results  returned  by  the  two  methods,  the 
perturbation  method  was  chosen  for  use  in  the  material  parameter  extraction.  Due  to  the 
large  air  gaps  and  incredibly  complex  nature  of  the  media,  it  is  assumed  that  the  Muller 
method  as  implemented  in  [8]  does  not  adequately  account  for  all  of  the  modes  present  in 
the  region.  Additionally,  the  Muller  method  did  not  allow  for  the  consideration  of  more 
than  2-3  modes  at  any  time  due  to  the  very  high  attenuation  present  in  .  This 
assumption  is  simple  at  best  and  requires  greater  exploration  to  confirm. 
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Table  4.1:  Comparison  of  the  wavenumbers  in  an  empty  eoaxial  line  and  a  line  partially 
filled  with  dieleetrie  with  a  123.5  mil  outer  air  gap.  The  frequeney  for  f  is  1  GHz.  The 
Muller  and  Perturbation  method  returned  nearly  identieal  results. 


Empty  Coax  Line 

Mode 

Y 

TEM 

0  +  21.0j 

TMoi 

771.0 

TMo2 

1552.1 

TMo3 

2331.2 

TMo4 

3109.7 

Perturbation  /  Muller’s  Method 

Mode 

7 

TMoo 

0.07  +  23.4) 

TMoi 

876.2-3.7) 

TMo2 

1606.7-2.7) 

TMo3 

2204.9  +  5.5) 

TMo4 

3020  +  2.7) 

Table  4.2:  Comparison  of  the  wavenumbers  returned  via  the  perturbation  and  Muller 
methods  of  a  line  eontaining  WaveX  with  a  139  mil  outer  air  gap.  The  frequeney  for  /  is 

1  GHz.  The  material  properties  used  for  these  ealeulations  are  =  83.29  +  749.98)  and 
(ju^  =  4.67  -  jl .11).  The  differenee  between  the  two  methods  indieate  that  Muller’s 
Method  perhaps  misses  several  modes  that  are  found  by  the  Perturbation  method. 


Perturbation  Method 

Muller’s  Method 

Mode 

7 

7 

TMoo 

11.4  +  32.1) 

11.4  +  32.1) 

TMoi 

882.7+  1.5) 

2,953.8  +  17.9) 

TMo2 

1,773  +  2.7) 

9,234.9  +  4.0) 

TMo3 

2,650.5  +  11.9) 

15,440.4-0.1) 

TMo4 

2,953.8  +  17.9) 

22,252.9  -  4.4) 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


The  focus  of  this  thesis  was  the  accurate  determination  of  the  complex 
permittivity  and  permeability  of  thin  samples  of  materials  which  partially  fill  a  coaxial 
transmission  line.  The  primary  method  used  to  determine  these  properties  was  the  modal 
analysis  method  developed  in  [8].  This  method  takes  into  account  higher  order  modes 
excited  by  air  gap  discontinuities  and  should  cause  the  results  to  converge  to  the  actual 
material  parameters  as  increasing  numbers  of  modes  are  considered.  This  assumes  that 
the  axially  symmetric  higher  order  modes  are  the  only  significant  source  of  error  in  the 
measured  values.  Other  methods  such  as  perturbation  techniques  and  a  low  frequency 
stripline  were  briefly  investigated  as  well  for  comparison  purposes. 

Both  a  low  loss  dielectric  material  and  high  loss  magnetic  material  were 
investigated  in  this  thesis.  The  dielectric  was  used  to  prove  the  effectiveness  of  the 
modal  method  with  large  air  gaps.  The  intent  was  to  utilize  a  2D  complex  root  search 
using  two  independently  obtained  values  of  82^  to  solve  for  and  ju^  since  it  was 
unknown  whether  the  S'jj  magnitudes  would  be  large  to  return  accurate  results  due  to  the 
large  air  gap.  It  turns  out  that  that  the  values  for  the  dielectric,  although  very  small, 
were  good  enough  to  be  used  in  the  2D  root  search  to  solve  for  and  ju^ .  This  technique 
was  then  extended  to  measured  samples  of  a  20  mil  thick  sample  of  the  magnetic  material 
WaveX.  The  WaveX  material  lacks  well  quantified  permittivity  and  permeability  data, 
so  the  low  frequency  stripline  method  was  briefly  explored  to  provide  data  to  compare  to 
the  modal  method.  The  development  of  this  technique  is  shown  in  Chapter  3  along  with 
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the  modal  analysis  method  development.  The  results  of  both  methods  are  reported  in 
Chapter  4. 

5.1  Future  Research 

There  are  several  avenues  that  could  be  explored  to  further  the  research  contained 
within  this  thesis. 

5.1.1.  Air  Gap  Sensitivity  Analysis.  This  research  considered  both 

dielectric  and  magnetic  materials  containing  outer  air  gaps  that  accounted  for  over  75% 
and  87%  respectively  of  the  cross  section  of  the  line.  The  modal  method  performed  well 
for  the  dielectric  sample  considering  such  a  large  air  gap,  but  performed  very  poorly  for 
the  WaveX  material.  The  performance  of  the  dielectric  sample  was  expected  based  upon 
results  of  previous  research  by  Fehlen.  The  MagRAM  samples  considered  in  [8] 
contained  only  a  9  mil  outer  air  gap  representing  only  approximately  1.5%  of  the  cross 
section  of  the  line.  The  method  performed  well  for  this  scenario,  but  still  contained  some 
error,  especially  at  the  lower  frequencies.  It  is  not  surprising  then,  given  the  very  large 
air  gap  considered  in  this  research,  that  the  method  would  perform  poorly  for  the  WaveX 
samples. 

One  avenue  of  future  exploration  could  be  the  sensitivity  of  the  modal  method  to 
the  size  of  the  air  gap  for  lossy  magnetic  materials.  The  use  of  a  more  well  quantified 
material  is  recommended  as  one  of  the  challenges  associated  with  WaveX  was  the 
availability  of  independently  determined  values  of  and  . 
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5.1.2.  Consideration  for  Anisotropic  Material.  The  modal  method 

development  performed  in  this  thesis  assumes  that  the  MUT  is  simple  media  (i.e.  linear, 
homogenous,  and  isotropic).  As  it  turns  out,  the  WaveX  samples  are  highly  anisotropic; 
therefore,  some  of  the  assumptions  made  in  the  field  developments  could  be  invalidated. 
A  closer  examination  of  the  effects  on  the  field  developments  for  an  anisotropic  material 
would  be  beneficial  beginning  with  Maxwell’s  equations. 
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APPENDIX  A:  MAXWELL  X  EQUA  LIONS 


The  comprehension  of  electromagnetic  wave  theory  and  the  mode-matching 
technique  begins  with  the  understanding  of  Maxwell’s  equations  which  are  summarized 
here  based  on  the  developments  from  [2,  6,  10].  This  understanding  is  the  basis  upon 
which  all  field  developments  for  waveguides  are  built. 


VxE-  -M  -  jcoB 

(A-1) 

YxH  =  J  +  jaD 

(A-2) 

(A-3) 

'^■B  =  q,nv 

(A-4) 

The  two  curl  equations  known  respectively  as  Faraday’s  Law  and  Ampere’s  Law  relate 
the  electric  (E)  and  magnetic  (H)  field  intensities  to  the  corresponding  magnetic  current 
and  flux  densities  (M  ,  B)  or  electric  current  and  flux  densities  ( J , .D )  as  a  function  of 
angular  frequency.  The  two  divergence  equations  known  as  Gauss’s  Law  for  Electric  and 
Magnetic  fields  relate  the  electric  or  magnetic  flux  densities  to  the  electric  or  magnetic 
charge  densities.  This  research  assumes  that  all  materials  of  interest  are  linear, 
homogeneous,  and  isotropic;  therefore, 

D  =  €E  (A-5) 

B  =  /jH  (A-6) 

with  £•  and  //  being  the  permittivity  and  permeability  of  the  material.  Making  this 
substitution  into  Maxwell’s  Equations  presents  them  in  their  simplified  form  as 
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(A-7) 


VxE-  -M  -  jcoiiiH 
V  X  //  =  J  +  jojsE  (A-8) 

W-E  =  ^  (A-9) 

£ 

V-H  =  ^.  (A-10) 

M 

Maxwell’s  equations,  (A-7)  through  (A-10)  are  first-order,  eoupled,  partial  differential 
equations.  They  eannot  be  solved  direetly  to  obtain  the  electrie  and  magnetic  fields; 
however,  they  can  be  used  to  form  second-order,  uncoupled,  partial  differential  equations 
which  are  readily  solvable  for  the  fields  in  the  system. 


67 


APPENDIX  B:  HELMHOLTZ  WAVE  EQUATIONS  &  VECTOR  POTENTIALS 


The  Helmholtz  wave  equation  for  the  electric  fields  is  formed  by  taking  the  curl 
of  both  sides  of  Faraday’s  Law  (A-7) 

Vx-VxE  --VxM  -jcojuVxH.  (B-1) 

Using  the  vector  identity 

VxVxH  =  v(V-h)-V^H,  (B-2) 

and  Ampere’s  Law  (A-8),  equation  (B-1)  now  becomes 

y{y-E)-V^E  =  -VxM-Jco^(J  +  jcosE).  (B-3) 

Applying  Gauss’s  Law  (A-9)  and  simplifying,  Helmholtz’s  electric  field  wave  equation 
becomes 

V^E  +  k^E  -  VxM  +  jcofjJ  +  —  (B-4) 

s 

where 

=  co^sii.  (B-5) 

Ampere’s  Law  can  be  treated  in  the  same  manner  yielding  Helmholtz’s  magnetic 
field  wave  equation 

V^H  +  k^H^-VxJ  +  jo)sM  +  -q^^.  (B-6) 

A 

Assuming  no  electric  or  magnetic  sources  are  present  in  the  region,  (B-4)  and  (B-6) 
reduce  to  the  source-free  electric  and  magnetic  field  Helmholtz  wave  equations: 

V^E  +  k'^E^Q  (B-7) 
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V^H  +  k'^H^Q. 


(B-8) 


The  fields  ean  now  be  solved  for  using  vector  potentials  and  the  separation  of 
variables  technique  described  in  [2],  The  wave  equation  of  the  form  of  (B-7)  or  (B-8) 
can  be  expressed  in  terms  of  cylindrical  coordinates  for  a  coaxial  system: 

VV(a^,2)  +  ^V(p,^,^)  =  0.  (B-9) 


When  expanding  (B-9)  in  cylindrical  coordinates,  the  wave  equation  becomes 


lA 

pdp 


V  ^P) 


+ 


1  ay 

d<i>^  ^ 


+  k^y/  -  0 


(B-10) 


where  ^  is  a  scalar  function  that  can  represent  a  field  or  vector  potential.  Assuming  that 


\f/{^p,(p,z^  has  a  separable  solution  of  the  form 

y/{p,(p,z)  =  f(p)g{(p)h(z), 

equation  (B-10)  can  be  rewritten  as 

1  f  I  I  df  11  d^g  1  d^h 
f  dp^  f  p  dp  g  p^  dcj)^  h  dz^ 


(B-11) 


(B-12) 


Since  the  next  to  last  and  last  term  on  the  left  side  of  (B-12)  are  functions  of  a  single 
variable,  they  must  be  constants  and  are  determined  to  be 


1  d^g  2 

- ir  =  —m 

g  df 


(B-13) 


1  d^h  ,  2 
— r  =  -k; 
h  dz^ 


Substituting  these  results  into  (B-12)  yields  the  reduced  equation 


A  A/ 

/  dp^ 


+ 


P  df 
f  dp 


+  +  p^  {k^ 


(B-14) 


(B-15) 
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Rearranging  the  terms  of  (B-15)  and  defining 


ki=e-ki 


(B-16) 


results  in  the  following  equations  defining  f,  g,  and  h  respeetively: 


d^g 

d(lP- 


dz^ 


■■-m  g 


=  -kth 


(B-17) 


(B-18) 


(B-19) 


Equation  (B-17)  is  readily  identified  as  the  Bessel  differential  equation  which 
means  the  general  solution  will  be  a  linear  combination  of  Bessel  or  Hankel  functions. 
For  the  coaxial  line  in  this  research,  the  field  will  be  a  standing  wave  in  the  p  direction 


yielding  the  solution 


f{p)  =  AJ^{kn)  +  BY^{kn) 


(B-20) 


where  and  represent  Bessel  functions  of  the  first  and  second  kind  respectively  of 
order  m. 

The  solutions  to  (B-18)  and  (B-19)  are  well  documented.  For  the  case  of  the  z  directed 
coaxial  waveguide,  the  solutions  to  these  are 

g(^)  =  Ccos(m^)-l-Z)sin(m^)  (B-21) 

h{z)  =  Ee-^^  +Fe^^  (B-22) 

where 

Y  =  jK-  (B-23) 
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The  solution  to  the  Helmholtz  wave  equation  (B-9)  can  now  be  expressed  by  a 


product  of  the  individual  solutions  of  (B-20),  (B-21),  and  (B-22) 

¥m,n  iP^  =  \_KJm  i^pp)  +  {k^p  „  COS  „  sin  (m^2l)]  x 

(B-24) 

where  subscripts  m  and  n  are  added  for  denoting  unique  solutions  and 


(B-25) 
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APPENDIX  C:  FIELDS  IN  A  COAXIAL  WAVEGUIDE 


C.l  TM  fields  in  a  Coaxial  Guide 

For  a  linear,  isotropic,  homogenous  material,  the  electric  and  magnetic  vector 
potentials  (F  and  A  respectively)  are 


F  =  0 


(C.l) 


The  magnetic  potential  A  takes  the  form  of  the  potential  function  (B-24)  and  becomes 


4  =  [4,,,4  (y  P)  +  4,.,  4  (y  p)]  [c„,„  cos  (m!«)  +  D,,, ,  sin  (m«»)]  X 

+  4,,„c’"]  (C-2) 

The  potentials  from  (C.l)  are  then  substituted  into  the  following  field  relationships: 


E  =  ^^\eA  +  vN-A]\- 


VxF 


JCOSJU 


+  -^rk^F  +  v(V-F)1. 


(C-3) 


//  jas/d 

Solving  (C-3)  using  the  magnetic  vector  potential  from  (C-2)  yields  the  generating 
functions  for  the  TM^  fields  present  in  the  coaxial  waveguide: 


jcosfd  dpdz 


E,=- 


1  1 
jcos/u  p  dpdz 


H  =11^ 
^  p  p  d^ 


p  dp 


E=- 


1 


f  a2 


jcosp 


dz^ 


-+r 


A^ 


i/  -0 


Filling  in  the  magnetic  vector  potential  into  these  generating  functions  yields  the  TM^ 
field  components  in  the  guide: 
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-koT 


jcos/d 


E  = 


Aj'n,  ikpp)  +  BY'^  {kpP)\  [c  cos  (m^^)  +  Z)  sin 
=  ■  \_AJ„, {k  p)  +  BY^ {k  p)\ [C sm{m(l))-D cos 

JCOS/up^  ^  JL  V  /J 

\_AJm  ikpP)  +  BY^  {kpP)\  [c  cos  (m^)  +  Z)  sin 
ikpP)  +  BY^  {kpp)\ [C  sin  (m^)  -  D  cos 
^Jm  ikpP)  +  BY'^  {kpP)\ [C  cos  (m^^)  +  D  sin 


Ee-^^-Fj^ 


JCOSfJ. 

m 


PP 

-k^  r 


Ee-r^-Fj^ 

Ee-^^+Fj^ 
Ee-y^+Fe^^ 
\Ee-^^^Fj^ 


After  the  field  eomponents  are  known,  boundary  conditions  at  both  the  inner 
(p=a)  and  outer  (p=c)  walls  of  the  waveguide  must  be  applied.  For  this  scenario,  both 
walls  are  eonsidered  to  be  PEC;  therefore,  this  requires  that  the  tangential  eleetrie  fields 
at  both  walls  must  be  zero. 


^Jmikpa)  +  BY^{kpa)  =  ^ 
=>  AJ^{kpC)  +  BY^{kpC)  =  0 


A  =  -B 


^mikpO) 

Jmikpa)' 


Solving  for  the  constant  A  at  the 
inner  wall  yields 


This  is  then  substituted  into  the  field  equations  produeing 
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where 


mB, 


BJm  (kpnO) 
-I 


—  jFo(^A:^„yC>^[Csin(m^2i)-£)cos(m^2i)][^£'e  +Fe^’'^ 


= - y— (^^„p)[Ceos(m^2i)  +  £)sin(m^2i)] 


Ee~^’'"  + 


K  [kpP]  =  [Jm  {ctp)  Yo  {aa)-J,  {aa)Y^  {ap)\. 

Assuming  the  waveguide  is  axially  symmetric  (m=0)  and  only  accounting  for  the 
forward  propagating  fields,  the  field  components  reduce  to 
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with  the  Eij,  and  Hp  fields  going  to  zero  due  to  axial  symmetry.  The  characteristic 

equation  needed  to  solve  for  Yn  is  produced  by  applying  the  boundary  conditions  at  the 
outer  wall  and  yields 


(>^pnc)  Yo  (k^„a)-Jo  (k^^a)  (k^^c) 


=  0 


(3.4) 


where  +  Yn  .  For  the  TEM  mode,  Yn  -  ^  and  from  [13],  the  value  for  Yn  for 


modes  is  well  estimated  by 
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Yn  = 

where  c  and  a  represent  the  outer  and  inner  radii  of  the  waveguide  [13].  This  is  used  as 
the  initial  guess  when  numerically  solving(3.4)  for  the  actual  value  of  Yn  for  higher  order 
axially  symmetric  TM  modes. 

C.2  TEM  Fields  in  a  Coaxial  Guide 

The  dominant  field  present  in  a  fully  filled  coaxial  waveguide  is  the  TEM^  field. 
This  field  is  also  sometimes  referred  to  as  the  TMoo  mode  and  can  be  developed  from  the 
axially  symmetric  (m=0)  TM^  fields  which  were  previously  derived.  These  axially 
symmetric  modes  allow  for  the  possibility  of  =  0  which  is  found  by  examining  the 

E^{^p  =  a,b)  =  Q  boundary  condition.  This  is  referred  to  as  the  n=0  mode.  In  this  case, 
kz=k  since  kp=0. 

The  TEM  fields  are  developed  by  taking  the  limit  as  ^  0  of  the  TM  fields 

developed  above.  The  small  argument  approximations  and  logarithmic  identity  shown  in 
Table  C.l  are  applied  to  yield  the  resulting  TEM  fields. 


Table  C.l;  Small  Argument  and  Eogarithm  Identities 
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Summarily,  axially  symmetric  TM^  fields  in  a  coaxial  waveguide  reduce  to  the  dominant 


TEM  fields 


where 


TTCOSjil 


(C.5) 


(C.6) 
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APPENDIX D:  GTRI  WAVEX RESULTS 


Multiple  samples  of  WaveX  (WX-A020)  were  measured  by  researchers  at  GTRI. 
The  low  frequency  measurements  were  taken  using  an  impedance  analyzer  fixture  while 
the  high  frequency  measurements  were  taken  using  a  focused  beam  system.  Both  the  in¬ 
plane  and  out-of-plane  complex  permittivity  was  extracted  along  with  the  in-plane 
complex  permeability.  The  “Old”  samples  correspond  to  a  WaveX  sample  manufactured 
in  the  2007  timeline  while  the  “New”  samples  correspond  to  a  WaveX  sample 
manufactured  sometime  in  201 1 .  The  presumption  is  that  manufacturing  processes  have 
improved  between  the  two  samples  explaining  the  differences.  The  extracted  parameters 
are  displayed  in  Figures  D.l  -  D.3. 


Figure  D.l;  Extracted  out-of-plane  permittivity  of  “Old”  WaveX  sample 
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Figure  D.2;  Extracted  in-plane  permittivity  of  “Old”  and  “New”  WaveX  samples 
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Figure  D.3;  Extracted  in-plane  permeability  of  “Old”  and  “New”  WaveX  samples 
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